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Introduction 

The Tschauner-Hempel equations are widely used to correct the separation distance drifts between a 
pair of satellites within a constellation in highly elliptical orbits [1]. This set of equations was discretized in 
the true anomaly angle [1] to be used in a digital steady-state hierarchical controller [2]. This controller [2] 
performed the drift correction between a pair of satellites within the constellation. The objective of a 
discretized system is to develop a simple algorithm to be implemented in the computer onboard the 
satellite. The main advantage of the discrete systems is that the computational time can be reduced by 
selecting a suitable sampling interval. For this digital system, the amount of data will depend on the 
sampling interval in the true anomaly angle [3]. The purpose of this paper is to implement the discrete 
Tschauner-FIempel equations and the steady-state hierarchical controller in the computer onboard the 
satellite. This set of equations is expressed in the true anomaly angle in which a relation will be formulated 
between the time and the true anomaly angle domains. 

Discrete Tschauner-Hempel Equations and Control System 

In reference [2], the discrete TH equations for a perturbed motion that defines the motion of a pair of 
satellites in an elliptical orbit about the Earth can be described as, 

y(k + 1) = A(k)y(k) + B(k)u(k ) + y/ [y(£)] + r(/(fc)) 

where, 

y(k)=[y l (k) y 2 (k) y 3 (k ) y\(k) y 2 (k ) y^A:)]; ()' = rf (^/ ; >>,(*)=( 1 + ecos f)xj(k) 

x, is positive against the motion of the spacecraft, x 2 is positive along the radial direction, and x, is 
positive when the right handed system is completed. The y equations are used as a mathematical 

transformation. i/j(k) defines the perturbation due to the oblateness of the Earth as explained in Ref. [1]. 
r ; ( f(k )) is the solar pressure force defined by the Karmyov theory [1]. The maneuvering spacecraft is 

assumed to have an applied thrust vector along the reference coordinate system, and the reference 
spacecraft is initially assumed to be acted on by a Newtonian gravitational force directed toward the center 
of the Earth. In this work, the satellites are assumed to be opaque right circular cylinders where the incident 
solar rays can affect the motion of the satellites. 

The steady-state hierarchical control scheme was implemented with the cost function defined in Ref. 
[2]. In steady-state [2], the control problem can be reduced and defined in the following form, 
y{k + 1) = [ A{k) - B(k)K(k)~ B{k)K v [k%{k] + [a/7 - BK p (k) jy/ \y{k)\ + r(/(*))] 
where A(k)= I + A fA(k), and B(k)= AfB(k) where I is a 6 X 6 identity matrix. A is the state matrix with 
dimensions 6x6, B(k) is the control matrix with dimensions 6x3, and i// [}’(/c)] is a 6 x 1 matrix. K(k) is the 
control gain, K v (k) is the gain due to the nonlinear terms in the y/(y(k )) column vector, and K p (k) is the 
gain due to the sum of the perturbations (y/(y(£)) + r(/(£))). 
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Implementation of the Control Problem 

In reference [2], a block diagram for the steady-state hierarchical control scheme was developed. This 
block diagram can be easily modeled in the computer onboard the satellite. The only problem is to relate 
the time, in which the thrusters will be fired, to the true anomaly angle. This relationship can be obtained 
through the angular momentum equation used in reference [4]. This equation shows that the true anomaly 
angle is an increasing function of time that can be used to determine the time. By developing an interface 
between the time, the sample separation distance and velocity, and the block diagram of the steady-state 
controller, the computer onboard the satellite can perform all the necessary computations to correct the 
separation distance conditions in a short period of time. 

Results 

Figure I shows the control solution for the drift correction [2] for a set of initial conditions when the 
eccentricity is equal to 0.8. In 
Figure 1, the separation 
distance and velocity drifts are 
corrected before the pair of 
satellites reaches the perigee 
point. In addition, the control 
effort shows a minimum 
consumption of thrust. In this 
problem, the time scale has 
been obtained from the angular 
momentum equation. This time 
scale shows when the thrusters 
should be fired to correct the 
separation distance drifts. 

Conclusion 

The present work 
contributes for the first time 
the implementation of the 
Tschauner-Hempel equations 
and the stead y-state 
hierarchical control scheme in 
the time domain. This discretized system can be easily implemented in the computer onboard the satellite 
with a simple control algorithm. The computational process will be performed in a short period of time 
because of the relationship determined between the time and the true anomaly angle. In this way, the 
controller can perform the desired maneuvers for the correction of the separation distance and velocity 
drifts for a pair of satellites. 
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